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ABSTRACT 

We provide predictions for the J/tp coherent production in the peripheral heavy ion 
coUisions at LHC and RHIC using the leading twist model of nuclear shadowing based on 
the QCD factorization theorem for diffraction and the HERA hard diffraction data. We 
demonstrate that for LHC kinematics this model leads to a bump-shape distribution in 
rapidity which is suppressed overall as compared to the expectations of the color transparency 
regime by a factor ~ 6. This is a significantly larger suppression than that expected within 
the impact parameter eikonal model. Thus we show that the interaction of spatially small 
wave package for which the total cross section of interaction with nucleons is small is still 
strongly shadowed by nuclear medium in high energy processes. 

1 Introduction 



Interaction of small size color singlet objects with hadrons is one of the most actively studied 
issues in high-energy QCD. In perturbative QCD (similar to QED) the total cross section of 
the interaction of such systems with hadrons is proportional to the area occupied by color 
within projectile hadron leading to the expectation of a color transparency phenomenon 
for various hard processes with nuclei. In the case of incoherent processes cross sections are 
expected to be proportional to the number of nucleons in the nuclei while in the case of 
coherent processes the amplitude is proportional to number of nucleons times the nuclear 
form factor. Possibility to approximate projectile heavy quarkonium as colorless dipole of 
heavy quarks can be formally derived from QCD within the limit when mass of heavy quark 
mg — > oo but XBj = ^roq/v is fixed and not extremely small 0. In this kinematics the size 
of heavy quarkonium is sufficiently small to justify applicability of PQCD. 

Recently the color transparency (CT) phenomenon was observed at FNAL by E791 
experiment 0] which studied the coherent process of dissociation of a 500 GeV pion into 



1 



two jets off the nuclei. The measurement has confirmed a number of predictions of Q in- 
cluding the A-dependence, and the transverse and longitudinal momentum distributions of 
the jets. Previously the color transparency type behaviour of the cross section was observed 
also in the coherent J/ip photoproduction at (i?^) = 120GeV 0. 

A natural question is whether the color transparency will hold for arbitrary high ener- 
gies? Two phenomena are expected to work against CT at high energies leading to onset 
of a new regime which we refer to as the color opacity regime. One is the leading twist 
gluon shadowing. Indeed the QCD factorization for hard exclusive coherent processes with 
nuclei like 'y^ + A^ ^^Vector meson" + A implies that the cross sections are proportional 
to the square of the gluon parton density Q^) at small x which is screened in nuclei 

as compared to the nucleon: GA{x,Q'^)/AGNix,Q'^) < 1. This obviously should lead to 
a gradual disappearance of color transparency H. Another mechanism for violation of 
CT at high energies is the increase of the small dipole-nucleon cross section with energy 
oc Gn^XjQ"^). For sufficiently large energies this cross section becomes comparable to the 
meson-nucleon cross sections and hence one may expect a significant suppression of the hard 
exclusive diffractive processes like DIS diffractive production of vector meson and photo- 
production of heavy quarkonium states as compared to the CT scenario. However it seems 
that this phenomenon is beyond the kinematics achievable for the photoproduction of J/tp 
mesons at RHIC {x ^ 2 ■ 10~^, Q'^ ^ lOGeV'^) and probably even at LHC. 

It was suggested in p, |^ to look for color opacity phenomenon using J/ip (photo) elec- 
troproduction. This however requires energies much larger than those available at the fixed 
target facilities and would require use of electron-nucleus colliders. At the same time esti- 
mates of the counting rates performed within the framework of the FELIX study [§] have 
demonstrated that the effective photon luminosities generated in peripheral heavy ion colli- 
sions at LHC would lead to significant rates of coherent photoproduction of vector mesons 
including T in reaction 

A + A^ A + A + V. (1) 

As a result it would be possible to study at LHC photoproduction of vector mesons in Pb- 
Pb and Ca-Ca collisions at energies much higher than the range W^p < IT.SGeV covered 
at the fixed target experiment at FNAL [Q. Note that even current experiments at RHIC 
[W-yp ^ 25GeV) should also exceed this limit. 

Currently the theory of photoinduced processes in AA collisions is well developed, for 
the recent review see . Hence we can combine it with our previous studies of the coherent 
photo(electro)production of vector mesons to make predictions for production of J/ip in the 
process (|l]). Typical transverse momenta which are exchanged between two nuclei in the 
peripheral collisions which leave nuclei intact are much smaller than the typical transverse 
momenta in the coherent photoproduction of vector mesons. As a result for the cross section 
integrated over the momentum of the nucleus which emits the quasireal photons we can use 
the standard Weizsacker- Williams approximation. 

Hence the cross section of the vector meson production integrated over the transverse 
momenta of the nucleus which emitted a photon can be written in the convoluted form: 

""^^^ Zk''^^ = ^ / ^^T..9f-^'r-..^,*Aky (2) 
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Here k = 'jk^ is the photon momentum in the colhding frame(/c3 - momentum of photon in 
the rest system of emitting nucleus and 7 - Lorentz factor). 

For the quantity n{k, b) presenting the flux of photons with momentum k in the collider 
frame we used the simplest approximated formQ: 

n{k,b)) = ^^,X'[K!{X) + -K',{X)], (3) 

where Kq(X) and Ki{X) are modified Bessel functions with argument X = ^ and b is the 

impact parameter distance between centres of colliding nuclei. The factor TAA^b) accounts for 
inelastic interactions of the nuclei at impact parameters b < 2Ra- It can be approximately 
calculated as 

TAAib) = J d%,TA{bi)TA{b - 61), (4) 

00 

where TA{b) = J dz pA{z,b) is the usual profile function of the nucleus. In our calcu- 

—00 

lations we use the nuclear matter density pa(-2, b) obtained from the mean field Hartree- 
Fock-Skyrme(HFS) model, which describes many global properties of nuclei as well as the 
intermediate and high energy elastic proton-nucleus scattering and nucleus electromagnetic 
form factors. This indicates that the HFS model provides a good description of both proton 
and neutron distribution in nuclei and takes into account a small difference between the 
matter distribution and the charge distribution. 

The main subject of our interest in this paper is estimating the cross section of the process 
7^4 J/ipA 

- 1 dZ-^M^JMhH, (5) 

where 

s = AEf^k = AjkniN (6) 

is invariant energy for 'j — N scattering {En = 'yrriN is the energy per nucleon in the 
cm. of the nucleus-nucleus collisions), t = — is square of the vector meson transverse 
momentum. 



2 Coherent photoproduction of J/ip off nuclei 

Let us discuss the photoproduction amplitude 7 + ^4— > J/ip + A in more details. We are 
interested here in the W^p range which can be probed at RHIC and LHC In this situation 
interaction of cc which in the final state forms J/ ip is still rather far from the black body 
limit in which cross section can be calculated in the model independent way fl^. Several 
mechanism of coherent interaction with several nucleons were suggested for this process. We 
focus here on the the leading twist mechanism of shadowing. There exist qualitative 
difference between the mechanism of interaction of a small dipole with several nucleons 
and the case of a similar interaction of an ordinary hadron. Let us for example consider 
interaction with two nucleons. The leading twist contribution is described by the diagrams 
where two gluons are attached to the dipole. To ensure that nucleus remains intact in such a 
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Figure 1: Correspondence between diagrams for hard gluon induced diffraction off nucleon 
and shadowing for the vector meson production. 



process we need to attach colorless lines to both nucleons. These diagrams are closely related 
to the diagrams corresponding to the gluon diffractive parton densities which are measured 
at HERA (see Fig. |I|) and hence to the similar diagrams for the gluon nuclear shadowing 

As a result it was possible to express the quark and gluon nuclear shadowing for the inter- 
action with two nucleons in a model independent way through the corresponding diffractive 
parton densities using the Gribov theory of inelastic shadowing [|ll|] and the QCD factor- 



ization theorem for the hard diffraction [0. An important discovery of HERA is that 



hard diffraction is indeed dominated by the leading twist contribution and gluons play a 
very important role in the diffraction(this is loosely referred to as gluon dominance of the 
Pomeron). Analysis of the DESY diffractive data indicates that in the gluon induced pro- 
cesses probability of the diffraction is much larger than in the quark induced processes []^. 
The recent HI data on diffractive dijet production [|13] provide an additional confirmation 



of this observation. Large probability of diffraction in the gluon induced hard processes 
could be understood in the s-channel language as formation of color octet dipoles of rather 
large size which can elastically scatter with a rather large cross section. The strength of this 
interaction can be quantified using optical theorem and introducing 

= ^Q j-..ff{x,Q')/dt{t = 0) 
(Ttot{x,Q^) 

for the hard process of scattering of a virtual photon off the gluon field of the nucleon. An 
important feature of this mechanism of coherent interaction is that it is practically absent 
for X > 0.02 -h 0.03 and may rather quickly become important with decrease of x. The 
gluon virtuality scale which is relevant for the J/ip photoproduction is 3-4 GeV^ with a 



significant fraction of the amplitude due to smaller virtualities 0, |T4|. Hence we will take 
the gluon shadowing in the leading twist at = AGeV"^. Taking a smaller value of would 
result in even larger shadowing effect. We present numerical values of al^^iQ"^ = AGeV"^) 
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for two current models of the diffractive gluon densities(Fig. |^) which practically cover the 
range of parameterizations available in the literature. The HI parameterization leads to a 
more graduate onset of the contribution of the double interactions because in this model 
diffraction into masses with M^/Q^ < 1 is smaller. The dijet data of HI prefer this scenario 
though it seems that further measurements will be necessary to clarify the issue. So we 
keep both models for the further analysis. For a more detailed discussion of the current 



models of diffraction and of the resulting values of cxf rr see [ff5]. The effective cross section 



'eff 

aeff can be used to estimate relative importance of the interactions with > 3 0, which 
corresponds to account of diagrams of Fig. ^ in the quasieikonal approximation. As a 
result the t-dependence of the photoproduction turns out about the same as for the case of 
Glauber scattering of a projectile with cross section of interaction with a nucleon equal to 
(Jeff- The ratio of the photoproduction cross sections off nucleus and nucleon is expressed in 
the leading twist through the ratio of the skewed gluon parton densities. In the case of J/ ip 
production they are pretty close to the gluon density calculated at Xg// = (xi +X2)/2 where 
Xi,X2 are light cone fractions carried by exchanged gluons and x = Xi—X2. In the DIS limit, 
or for tuq — >■ oo, one finds [|l^ X2 -C x and hence Xg// ~ x/2. However in the case of J/ip 
photoproduction Fermi motion effects lead to X2/X ~ 0.3 — 0.5 and hence to Xg// ~ x. 

In principle the multiple eikonal type rescatterings (at fixed transverse separations) due 
to multiple gluon exchanges - see Fig. § (the impact parameter eikonal rescattering 
model) could also result in suppression of the vector meson production. Though validity 
of this approximation is hard to justify in QCD the model calculations suggest that this 
effect is not small numerically 0. However, it is still significantly smaller than the leading 
twist shadowing (at least for x < 0.001) which we find in our calculations. Note in passing 



that if one would consider the gluon shadowing using phenomenological models [|T^, 
where shadowing for gluons was assumed to be equal to that for quarks at low normalization 
scale one would find comparable suppressions due to the leading twist gluon shadowing 
and the eikonal rescatterings. This would make extraction of the gluon shadowing from 
measurements of photoproduction of J/ ijj highly problematic. 



3 Numerical results and discussion. 

Having build the quasieikonal model for the amplitudes of the interaction of cc pair with 
several nucleons we can now calculate the amplitude of scattering off nuclei. We have demon- 
strated in that the amplitude in this approximation has the same structure of the rescat- 
tering terms as the Glauber model with atot substituted by <Jeff- Hence we can use the 
optical limit of the Glauber model[|l8| to calculate the cross section of J/ip photoproduction 



da^A-^VAis,t) _ daj]\f^vNis,t = 0) 



dt 



dt 



d%dze''^'-^p{b,z)e''''' -e 



-\a,ff{Mlrs) ] p(b,z')dz' 



Here the exponential factor with qi = my/2k accounts for finite longitudinal distances in 
the transition 7 — > (finite longitudinal momentum transfer). The forward elementary cross 
section for photoproduction of J/ ip meson on nucleon was taken using the fit to experimental 



5 




6 



Figure 3: Leading twist diagrams for the production of J/ip off two and three nucleons. 




Figure 4: Typical diagrams for the higher twist eikonal interactions of a small dipole with 
two nucleons. 
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data presented in |T^ (this is preferable to using the theoretical calculations which for 
photoproduction of J /ip have theoretical uncertainty of the order of two). 
We focus here on the distributions over rapidity: 



.--ml^-ln^. (9) 



In Figs. ^, ^, and Fig. |^ we present the differential cross sections both including effects 
of gluon shadowing and without gluon shadowing (impulse approximation) for lead-lead 
peripheral collisions at LHC and gold-gold collisions at RHIC. 

One can see that on the top of the overall suppression of the cross section the gluon 
nuclear shadowing leads to a significant modification of the shape of the rapidity distribution 
as compared to the impulse approximation. Bumps near the edges of rapidity distribution 
are due to a sharp increase of the effective cross section found in the calculation based on 
the Alvero et. al. model of the gluon diffractive density in the region of Bjorken x close to 
X ^ 10~^ . A bump in the center region of y arises due to the drop of the photon flux as can 
be seen from Fig. ^. These effects are weakened for the HI parameterization which leads to a 
more gradual increase of aefj and they disappear in the impulse approximation. It is also of 
interest that in the LHC kinematics we are sensitive to the cross section of photoproduction 
at W^p up to a factor of three larger than W^p corresponding to production at y = 0. Hence 
the measurements will actually probe the J/ ip photoproduction at the energies beyond those 
reachable at HERA in electron-nucleus mode. 

In the RHIC kinematics we find even more nice picture in the case of gold-gold collision. 
The decrease of cross section as a function of rapidity due to the shadowing is combined 
with drop of the photon flux in the same region of rapidities. This results in a narrow dip 
at y = which is very sensitive to pattern of onset of the gluon shadowing. The test of this 



prediction will be feasible at RHIC since the rates of J/ip production are pretty high pT[| . 

At LHC the total cross section calculated with accounting for the gluon shadowing effects 
is a{PbPb J/i) + PhPh) 14 mb and the value a{PbPb J/ip + PhPh) 85 mb 
was obtained in the impulse approximation. Hence, we have found a strong suppression 
of the J/ip yield for this case as it was predicted on the base of rough estimates in[0. In 
kinematic of RHIC the effect of suppression due to the gluon shadowing is rather small 
and a{AuAu ^ J/ip + AuAu) ~ 0.320 mb while in the impulse approximation a{AuAu 
J/ip + AuAu) ^ 0.360 mb. 

Let us compare our results with two other calculations of the reaction (1). The first rather 



detailed calculation of the coherent process A A —^A + V + A has been reported in Ref.|21 
To evaluate nuclear shadowing effects in the total J / if) A cross section the vector dominance 
model, classical mechanics formulae (accounting for the elastic rescatterings of vector meson 



only) have been used in [^. On the contrary our calculation uses eikonal approximation 
where inelastic shadowing effects dominate. Really dgj^ derived from the diffractive gluon 
densities includes both the elastic and inelastic shadowing. It is also assumed in Ref. |^ that 
the t-dependence of the cross section is cx |FA(t)p (where FAif) is the nuclear form factor) 
while account for the rescattering effects (eq. ^ leads to somewhat steeper t-dependence. For 
the RHIC energies whenever comparison is possible the results of Ref. PT| are pretty close to 



ours. This is because nuclear shadowing effects are a small correction for the photoproduction 
of J /if) in the kinematics of RHIC. For the LHC kinematics we obtained cross section for the 
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Figure 5: The rapidity distribution for the J /if) production in lead-lead peripheral collisions 
at LHC Solid line - production by two-side beams, dashed - production by one-side beam 
only. Calculations were performed with cr///^ based on Alvero et.al. parameterization of 
the gluon diffractive density. 
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Figure 6: The rapidity distributions for the LHC lead-lead peripheral collision J /ip coherent 
production calculated with uif^^ based on HI (solid line) and Alvero et al (dashed line) 
parameterizations of gluon density and in the impulse approximation (dot-dashed line) 
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Figure 7: The same as in Fig. ^ but for gold-gold collisions at RHIC. 
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coherent J/ip production significantly below the value of [^| (a factor of two for lead -lead 



collisions). The difference is because approach used in significantly underestimates the 
strengths of multiple interaction of cc pair with the nucleus. Also rapidity distributions for 



lead-lead collisions for which we find an interesting shape were not considered in Ref.pT 



After this study was nearly completed a report has appeared |^| where it has been 
suggested to use the coherent J /ip production in the peripheral ion- ion collisions to measure 
shadowing of gluon densities in nuclei. The analysis in p2| is based on the factorization 
theorem of ^, the t-dependence of the coherent 7^4 J/ipA cross section has been 
approximated as F'^(t) and three sets of the gluon distributions has been used. Calculations 
with the GRV gluon distribution is effectively equivalent to the impulse approximation. Two 



others model accounts for the nuclear shadowing. The distribution of [|T7| assumes the same 
shadowing for gluons as for quarks which is in variance with diffractive data from HERA. 
The second model |]23[ attempts to account for nonlinear QCD evolution in the gluon density 
and in this case application of the factorization approximation is hardly justified. It would be 
reasonable to expect that at least for J/ip production in the kinematics at RHIC where the 
shadowing effects don't influence essentially the total cross sections our numerical results 
should be qualitatively similar. However we found gross differences - our absolute cross 



sections are considerably larger as compared to that in |^ both for the kinematics at RHIC 
and at LHC. 

To conclude, we have demonstrated that heavy ion collisions at RHIC are sensitive to 
the onset of gluon nuclear shadowing while the measurements at LHC will allow to establish 
disappearance of transparency of nuclear matter for spatially small cc dipole at high energies 
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